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I Introduction 


Recently, much attention has been paid to the rare flavor-changing neutral current decay B —> K*{—> Kir)p + pT. 
One of the reasons for this was the first measurement of form-factor-independent angular observables performed by 
the LHCb Collaboration 0,i- It has been claimed that there is a 3.7cr deviation from the Standard Model (SM) 
prediction for one of the angular observables. Much effort has been spent to explain this deviation by invoking 
the effects of new physics (NP) (for example, see Refs. and references therein). The main emphasis of the 
above-mentioned papers was on the search for the physical observables that have low sensitivity to the form factors. 

In addition to the NP effects, the uncertainties related to the presence of the intermediate scalar resonance Kq 
decaying into Kir have been intensively discussed in the literature [lCHlfll . A detailed analysis of the B —> Kj(—> 
Kir)p + p~ decay in the higher kaon resonance region was done in Ref. TO]. In many papers, the Breit-Wigner form 
for the Kir mass spectra was used. However, this assumption cannot be justified for the broad scalar resonances like 
the A^q( 800 ) meson. The improvement of the description was done in Ref. [TTJ by invoking the chiral perturbation 
theory for the Kir interaction. This issue was also generalized to B s —> Kirlv in Ref. [l2| . 

As is well-known, short-distance physics is under control in the description of the rare B decays, whereas the 
effects of long-distance physics described by the hadronic form factors lead to large uncertainties since they involve 
nonperturbative QCD. The calculation of the B —> K* transition form factors have been performed in many theoretical 
approaches and models. We must mention some of them: light-cone QCD sum rules TtI . QCD sum rules 0 , the 

lattice-constrained dispersion quark model [19 [_j the simple dipole parametrization [20 1, perturbative QCD at large 

recoil region [2lj| . the relativistic quark model 22|, and the Dyson-Schwinger equations in QCD [23j . 

The B s and D s to Aq( 1430) transition form factors were calculated in Ref. (27| within an approach based on QCD 
sum rules. The form factors for the B —> ATq( 1430) transition have been evaluated in the light-front quark model 
[28| . The form factors of rare B —t Kq (1430)£ + £ _ decay were calculated in Ref. [29| within three-point QCD sum 
rules. The B -A S transition form factors have been investigated in the light-cone sum rules approach [3Cj. The 
transition form factors of R(H s )-mesons decay into a scalar meson were studied in Ref. [31] within the perturbative 
QCD approach. With these form factors, the decay width and branching ratios of the semileptonic B —> S£D and 
rare B -A SI + I~ decays have been calculated. The rare semileptonic decays B s —»■ [/ o (980), Aq( 1430)]£ + £ _ and 
B s — > [jo (980), Kq (1430)] vv were investigated in Ref. [H] in the framework of the three-point QCD sum rules. The 
B s — > /o(980) transition form factors were computed in Ref. 0 by using light-cone QCD sum rules at leading order 
in the strong coupling constant and an estimate of next-to-leading-order corrections. A QCD light-cone sum rule was 
also used to evaluate the B s —► S form factors and B s —t Slvt, M [I = e,/x,r) branching ratios in Ref. (H). The 
twist-3 light-cone distribution amplitudes (LCDAs) of the scalar mesons were investigated in Ref. [35j within the QCD 
sum rules. As an application of those twist-3 LCDAs, the B —► S transition form factors were studied by introducing 
proper chiral currents into the correlator. 

Recently, the B — S form factors for two scalar nonet mesons below and above 1 GeV were calculated in Ref. [36] 
by taking into account the perturbative 0(a s ) corrections to the twist-2 terms using the light-cone QCD sum rules. 
They were used in Ref. |37j to study the semileptonic B —i Sfvt and rare B —t S £ + £~, Svv decays. 

In the wake of exploring uncertainty in the full angular distribution of the B —>■ K*(—> Kir)n + pT decay caused 
by the presence of the intermediate scalar Kq meson, we perform the straightforward calculation of the B{B S ) —> S 
(S is a scalar meson) transition form factors in the full kinematical region within the covariant quark model. We 
restrict ourselves to the scalar mesons below 1 GeV: ao(980), jo(500), jo(980), and Ao(800) [241 ]. Actually, the 
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internal structure of these mesons is not yet well established (see Refs. [25), [26J for a review). We will use the simple 
qq interpretation of the low-lying scalar mesons in our calculation. The calculated form factors are used to evaluate 
the branching fractions of the decay B(B S ) S££, where £ = e, //, r. We compare our results with those obtained in 
other approaches. 

The paper is organized in the following manner. In Sec. II we give the necessary theoretical framework which includes 
the effective Hamiltonian, its matrix element between the initial and final states, the definition of the hadronic form 
factors and the helicity amplitudes. In Sec. Ill we briefly discuss our covariant quark model and calculate the form 
factors of the transitions B —> Sivi and B —> S£ + £~. Finally, we present our numerical results for the differential 
decay distributions and branching ratios. We compare our findings with the results of other approaches. 


II Effective Hamiltonian and form factors 

We start with the on-shell decays Bd —> (K, Kq , K*)£ + £~ which can be described by using the effective Hamiltonian 
for the b —> s transition [H), [39J. The effective Hamiltonian leads to the free-quark b —»■ sl + l~ decay amplitude: 

m (5 — >■ s£ + £~) = {cf {son) (e^e) + c 10 (so^b) (I w e) 

- (s i^ v (1 + 7 5 ) O) M } (1) 

where O^ = 7^(1 — q 5 ) is the weak Dirac matrix, A* = | V f { s Vtb | is the product of the Cabibbo-Kobayashi-Maskawa 
elements, and Cf = CV — Cs/3 — Cq. The Wilson coefficient Cf effectively takes into account (i) the contributions 
from the four-quark operators and (ii) the nonperturbative effects coming from the cc-resonance contributions which 
are as usual parametrized by a Breit-Wigner ansatz [40]: 
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( 2 ) 


where /1 is a scale parameter and m-i = ms- In what follows, we will not include the long-distance contributions 
coming from the J/ij} and i/j(2S) resonances [4C| and charm-loop effects [4l| . 

We specify our choice of the momenta as p\ = P 2 + k\ + k 2 with p\ = mf, p 2 = m 2 and kf = k 2 = m 2 where fci 
and k 2 are the £ + and £~ momenta, and mi, m 2 , ran are the masses of the initial meson H\ 1 the final meson H 2l and 
the lepton £, respectively. The matrix elements of the exclusive transitions B —> K(Kq)££ are defined by 


M(Hi H 2 ££) 


{Cf < HzlsOniH! > C fll £ 

Cio < H 2 \s0^b\H 1 > £^^q£ 

Cf < H 2 I sia^ (1 + 7 5 ) q v b\H 1 > *y/} 


( 3 ) 









where H x = B, H 2 = K(I<*). 

We define dimensionless form factors by 
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<H 2 (p 2 )\sO»b\H 1 (p 1 )> = Ffq 2 )P» +Ffq 2 )qd 

< H 2 (p 2 )\sia IJ ’ v q u (l+'Y 5 )b\H 1 (p 1 ) > = --j-- (P^q 2 - q^ Pq) F T (q 2 ) , (4) 

m i + m 2 

where P = Pi + p 2 and q = pi — p 2 . The matrix element in Eq Q is written as 

M (H, ^H 2 +U) = ^-^ {Tf flpl) + {h^d)} 

where the quantities Tf 1, are expressed through the form factors and the Wilson coefficients in the case of the spinless 
particle H 2 as 


T? = Ffp^+Ffqf 1 (i = 1,2), 

F[ 1] = Cf F + + Cf F t 2mb , 

m i + m 2 

pf = Cf F_ - Cf F t 2mb ^ , 

mi + m 2 q z 

J-£ 2) = C 10 F±. 

Respectively, the helicity form factors H l m are defined in terms of the invariant form factors as (44 

1 


m = 




(- PqF + + q 2 r _) 


H l ± =0 


H i = 2w ljP2| p 




( 5 ) 

( 6 ) 


The differential ( q 2 ,cos9 ) two-fold decay distribution may be written in terms of the bilinear combinations of the 
helicity amplitudes (see Ref. (42[). However, it is common in the modern literature to use the transversality amplitudes 
and A t defined in Ref. 0]- They are related to our helicity amplitudes by 
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- H { f) 
+ Hf) 


where the overall factor is given by 


T G% ( aAA 2 |p 2 |g 2 ^ l s 
.4 (27r) 3 \2n ) 12m 2 . 


( 7 ) 


where |p 2 1 = A 1 / 2 (mf, m|, q 2 )/2m\ is the momentum of the outgoing meson H 2 and v = fl — 4 m 2 /q 2 is the lepton 
velocity, both of which are given in the rest frame of the parent meson Hi. 
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The differential decay distribution then reads 

dT{Hi -> H 2 U) _ G 2 f faX t - 2 

dq 2 d {cos 9) 
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Integrating over cos 0 one obtains 

dT(Ui —> H 2 il) 
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( 8 ) 


(9) 


where we have introduced a flip parameter 8 u = 2 m 2 /g 2 . 

We also calculate the differential rates for the semileptonic Hi —> H 2 (.i/£ mode and rare Hi —> H 2 i/D decay. One 
has 
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The form factors F± are defined by Eq. whereas the function D u {xt) is given by 
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( 11 ) 


III The B — S transition form factors in the covariant quark model 

We calculate the B — S transition form factors in the covariant quark model. We briefly recall the basic features of 
this approach, which was formulated in its modern form in Ref. [44j] by taking into account the infrared confinement 
of quarks. 

The model is based on an effective interaction Lagrangian describing the coupling of hadrons to their constituent 
quarks. For instance, the coupling of a meson M{qiq 2 ) to its constituent quarks qi and q 2 is described by the nonlocal 
Lagrangian 






( 12 ) 


Here, Vm is the Dirac matrix, which is chosen appropriately to describe the spin quantum numbers of the meson field 
M{x). The vertex function Fm{x,x i,x 2 ) characterizes the finite size of the meson. To satisfy translational invariance 
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the vertex function has to obey the identity Fm{x + a,X\ + a,x 2 + a) = Fm(x,x i,x 2 ) for any given four-vector a. 
We use a specific form for the vertex function which satisfies the above translation invariance relation. One has 


2 

F m {x,xi,x 2 ) = d^jx -^WjXj) <1 >m 

i =1 


^(*1 - x 2 ) 


(13) 


where $>m is the correlation function of the two constituent quarks with masses m\ and m 2 . The variable Wi is defined 
by Wi = mi/{rrii + m 2 ), so that w± + w 2 = 1. We choose a simple Gaussian form for the vertex function < 1 > m(— k 2 ). 
The minus sign in the argument of < J>m(— fc 2 ) is chosen to emphasize that we are working in Minkowski space. One 
has 


$m(-^ 2 ) = exp(fc 2 /A^) 


(14) 


where the parameter A m characterizes the size of the meson. Since k 2 turns into —k^ in Euclidean space the form (1141) 
has the appropriate falloff behavior in the Euclidean region. We stress that any choice for is appropriate as long 
as it falls off sufficiently fast in the ultraviolet region of Euclidean space in order to render the Feynman diagrams 
ultraviolet finite. 

In the evaluation of the quark-loop diagrams we use the free local fermion propagator for the constituent quark, 


S q (k) 


1 

m q — fi — ie 


(15) 


with an effective constituent quark mass m q . 

The coupling constant gM in Eq. (1121) is determined by the so-called compositeness condition suggested by Wein¬ 
berg [45| and Salam j46| (for a review, see Ref. 0) and extensively used in our studies (for details, see Ref. 0). 
The compositeness condition requires that the renormalization constant Zm of the elementary meson field M(x) is 
set to zero, i.e., 


Zm — 1 — 5 , Mn , M (m| f ) — 0 


where W M (p 2 ) is the derivative of the mass operator corresponding to the self-energy diagram in Fig. [T] 


(16) 


qi 



M 


Q2 


FIG. 1: Diagram describing the meson mass operator. 


1/2 

To clarify the physical meaning of the compositeness condition, we recall that the renormalization constant can 
also be interpreted as the matrix element between the physical state and the corresponding bare state. For Zm = 0 it 
then follows that the physical state does not contain the bare one and it is therefore described as a bound state. The 
interaction Lagrangian (1121) and the corresponding free Lagrangian describe both the constituents (quarks) and the 
physical particles (hadrons), which are bound states of the constituents. As a result of the interaction, the physical 
particle is dressed, i.e., its mass and wave function have to be renormalized. The condition Zm = 0 also effectively 
excludes the constituent degrees of freedom from the space of physical states and thereby guarantees that there will 
be no double counting. The constituents exist in virtual states only. 

The covariant quark model was applied to evaluate the form factors of the B{B S ) —> P(V)— transitions in the full 
kinematical region of momentum transfer squared [49, 50| This approach was extended to describe the baryons as 
three-quark states (5l[ and the exotic meson X(3872) as a tetraquark [53 ] . 
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A similar approach based on the compositeness condition Z = 0 was recently developed in Ref. [53}. 

In this paper we evaluate the B—S transition form factors assuming that the scalar mesons below 1 GeV are ordinary 
two-quark states. Some remarks should be made before performing the calculations. The internal structure of the light 
scalar mesons is not yet well established (for review, see Refs. [25|, [26|] ). Since they have large decay widths it is difficult 
to distinguish them from background. There are interpretations of these objects as four-quark states and/or gluballs. 
Here, we describe the scalar mesons as two-quark states and evaluate the B — S form factors within our approach, but 
when we use the calculated form factors in the matrix element of the cascade decay B —> Kq(—> Ktt)£ + £~ we take 
into account the line shape of the Kq, which reflects the broad width of this resonance. One can also describe the 
scalar mesons as four-quark states in our approach, similar to the exotic meson X(3872) [52}; however, this is beyond 
the scope of this work. 

The SU(3) nonet of scalar mesons below 1 GeV can be written in the matrix form 


S = -= V S l X l 


X° = \ - I. 


i=0 


The physical scalar fields are related to the Cartesian basis in the following manner: 


S ± = j 5 (S 1 TiS 2 ), 

S+=-L(S* - iS 5 ), 

S7 = ^(S* + iS 5 ), 

S' = S° cos 9 S + s s sin 9 S . 


S° = S 3 , 

5 ° = 7f( 56 -* 57 )’ 

5° = ^(S 6 +*S 7 ), 

S = —S° sin 9s + S s cos 9s , 


(17) 


(18) 


where 9s is the octet-singlet mixing angle. The qSq vertex is then written as 
C-Sqq = qSq 

= S + ud + S~ du + S° (uu — dd) + S p us + S® ds + S~ su + 5° sd 

+ S'(^cosSs^(uu + dd) — sinds ss'j — S (^smds^(uu + dd) + cosSs ss'j , (19) 

where 5s = 9 — 9j, with the ideal mixing angle 9j = arctan (l/y/2). We will use the notation from Ref. [24} for the 
scalar mesons below 1 GeV: 

. = K* 0 { 800), I ( J p ) = A(0+), m K . {ao0 ) = 682 ± 29 MeV; 

. S' = / o (500), I G ( J PC ) = 0+(0++), m/ 0( 5 oo) = 400 - 550 MeV; 

. 5 = /o(980), I G (J PC ) = 0+(0++), to / o(980) = 990 ± 20 MeV; 

. S±- 0 = oj-°(980), I G (J PC ) = l-(0++), TO ao(980) = 980 ± 20 MeV. 

Moreover, we assume that Ss = 0, i.e., toj 0 ( 980 ) to ensure a pure ss state. 

The coupling constant gs in Eq. 113 is determined by Eq. msi), where is the derivative of the scalar meson 
mass operator, 


n s(p 2 ) = tr 


J_ f d 4 k 

2 p 2 ^ dp a J 4 n 2 i 
1 rd A k 

2 p 


I 4u 2 ! 


M a ) { 


2w ” Si(k + wip) S 2 {k - w 2 p) 
wi tr Si(k + wip) t/> Si(k + wip) S 2 {k — w 2 p) 


— w 2 tr 


Si(k + wip) S 2 (k — w 2 p) $S 2 (k — w 2 p) 


( 20 ) 
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By using the calculation technique outlined in Ref. |44|, 
analytic result for equal quark masses (m gi = m q2 = m q ): 


one can easily perform the loop integration. We give the 


IW) = 


1/A 2 1 

rah— 

0 0 


X 4{p 2 ^ 4 ™g + — [20 + t(l — 2a) 2 (1277i^ — p 2 )] - -4(1 - 2a) 2 (12 +p 2 f) + -^-(1 - 2a) Vj (21) 

v. 0'S CL a Clo * 


st 


z 0 = am q - a(l - a)p 2 , z x =(l - 2a) 2 p 2 

2 a s 


as = 2s + t , s = 


A 2 ' 


Note that in the case of A —> 0 the branching point appears at p 2 = 4to^. At this point the integral over t becomes 
divergent as t —> oo because zo = 0 at a = 1/2. By introducing an infrared cutoff on the upper limit of the scale of 
integration, one can avoid the appearance of the threshold singularity. 

Herein our primary subjects are the B — S transition matrix elements, which can be expressed via the dimensionless 
form factors defined in Refs. fTiil. [5(1 . The diagram corresponding to these matrix elements is shown in Fig. [2] 


B(pi) 

► 



<7i k + pi 


<73 
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= 7 e(i _ 7 5 ) i 1 + 7 5 ) 



- (fc + w 13 Pi) 2 ) 


- (fe + W23P2) 2 ) 


FIG. 2: Diagrammatic representation of the matrix elements in Eqs. (1221) and (1231) . 








One has 


(>%392](P2) I 92 <71 |B [9 - i9s] (pi)> 

iV c 3 b 3s J ( “ ( fc + Wl 3Pi) 2 ) ^ ( k + w 23 p 2 ) 2 ^ 


x tr 


S 2 (k+p 2 )O^S 1 (k+p 1 ) 7 5 S 3 (k) 


= F Bb (q 2 )P» + F Bb (q 2 )q^ , 


( 22 ) 


( 5 [ 93 9 2 ](P 2 ) I 92 (ia^q^l + j 5 )) q 1 \B [qiq3] ( Pl )) 

N c 3b 3s J ^TT^i ( ~ (^ + Wl3Pl) 2 ^ $S ( ~ (fc + W23P 2 ) 2 ^ 


x tr 


52 (fc + p 2 )icr^q v (\ + 7 s ) Si(k + pi) 7 s S 3 (fc) 
1 


mi + m2 


(q 2 P»-q-Pq») F BS {q 2 ). 


(23) 


Here, pf = m f, q± = b, q 2 = u, s, d , and 33 = s, d. Since there are three sorts of quarks involved in these processes, we 
introduce the notation with two subscripts, Wij = m qj /(m qi + m qj ) ( i,j = 1,2,3) so that Wq + Wji = 1. 

The first fit of the model parameters was done in the original paper [44[, where the infrared quark confinement 
was implemented for the first time. The leptonic decay constants (which are known either from experiments or from 
lattice simulations) have been chosen as the input quantities to adjust the model parameters. A given meson in 
the interaction Lagrangian is characterized by the coupling constant 3 b', the size parameter A# and two of the four 
constituent quark masses, m q ( m u = mj, m s , m c , mb). Moreover, there is the infrared confinement parameter A, 
which is universal for all hadrons. Note that the physical values for the hadron masses have been used in the fit. 
Therefore, one has 2 nn + 5 adjustable parameters for nn numbers of mesons. The compositeness condition provides 
nn constraints and allows one to express all coupling constants gn via other model parameters. The remaining 
nn + 5 parameters are determined by a fit to experimental data. The values of leptonic decay constants and some 
electromagnetic decay widths have been chosen as the input data. Several updated fits were done in Refs. [49l. [ 5 ( 1 ] ■ 
In this paper we will use the latest fit done in Ref. (54|. The fitted values of the constituent quark masses m q , the 
infrared cut-off A, and the size parameters A H are given by Eq. (EH) and Table Q] 


'Wlu/d 

m s 

m c 

m b 

A 


0.241 

0.428 

1.67 

5.05 

0.181 

GeV 



TABLE I: The fitted values of the size 

parameters A h 

in GeV. 


7T 

K D 

D s 

B 

Bs 

B c 

r/c 

Vb 


0.87 

1.02 1.71 

1.81 

1.96 

2.05 

2.50 

2.06 

2.95 


P 

Ul <j> 

J/V> 

K* 

D* 

D* 

B* 

B* s 

T 

0.61 

0.50 0.91 

1.93 

0.75 

1.51 

1.71 

1.76 

1.71 

2.96 


(24) 


Our form factors are represented as three-fold integrals which are calculated by using NAG routines. The results 
of our numerical calculations are well approximated by the parametrization 


F(q 2 ) = 


F( 0) 


1 — as + bs 2 ’ 


(25) 


We consider the following weak transitions: b — u (charged current), and b — d and b — s (flavor-changing neutral 
currents). The values of F( 0), a, and b are listed in Table HT1 
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1 

kC5 

to 

B-S 

A s = 0.8 

As = 1.5 

F+( 0) 

a+ 

b+ 

F+( 0) 

a+ 

b+ 

b — u 

B° d -a+( 980) 

0.144 

1.624 

0.585 

0.192 

1.433 

0.381 

b — u 

Bg — Kg + (800) 

0.138 

1.667 

0.674 

0.274 

1.258 

0.292 

b — s 

Bg — /o(980) 

0.141 

1.663 

0.651 

0.254 

1.269 

0.262 

b — s 

B° d -K5°( 800) 

0.191 

1.348 

0.407 

0.306 

0.988 

0.108 

b — d 

B d — /o(500) 

0.120 

1.448 

0.485 

0.210 

1.067 

0.155 


qi - <72 

B-S 

As = 0.8 

As = 1.5 

~F-( 0) 

a~ 

b_ 

~F-( 0) 

a- 

b _ 

b — u 

52 - 4 ( 980 ) 

0.049 

2.144 

1.196 

0.089 

1.723 

0.688 

b — u 

B°g-K*+{ 800 ) 

0.138 

1.727 

0.734 

0.268 

1.291 

0.310 

b — s 

to 

Co o 

1 

00 

o 

0.140 

1.761 

0.755 

0.253 

1.320 

0.295 

b — s 

B° d -KS°( 800) 

0.199 

1.406 

0.457 

0.296 

1.032 

0.129 

b — d 

B°d — /o(500) 

0.116 

1.504 

0.536 

0.191 

1.110 

0.180 


<7i - <72 

B-S 

As = 0.8 

As = 1.5 

F t { 0) 

CLT 


F t { 0) 

dT 

6t 

b — s 

B° - /o(980) 

0.165 

1.680 

0.667 

0.285 

1.276 

0.257 

b — s 

o 

o 

00 

o 
* o 

i 

o-s 

0) 

0.206 

1.367 

0.423 

0.306 

1.005 

0.113 

b-d 

B d — /o(500) 

0.124 

1.460 

0.496 

0.203 

1.080 

0.159 


TABLE II: The parameters of the fitted transition form factors F(q 2 ) at As = 0.8 GeV and As = 1.5 GeV. 


IV Numerical results and discussion 

We use the following set of SM parameters: Gp = 1.16637 x 10 -5 GeV -2 , rh c = 1.27 GeV, ra/, = 4.19 GeV, 
fht = 173.8 GeV, mw = 80.41 GeV, sin 2 9w = 0.2233, A * = |Vt s V' t [,| = 0.041, \V u b\ = 0.00413, a scale parameter 
At = rh b , and the Wilson coefficients Ci = -0.248, C 2 = 1.107, C 3 =0.011, C 4 = -0.026, C 5 = 0.007, C 6 = -0.031, 
Cy B == —0.313, Cg = 4.344, and Cio = —4.669. We take the average values of the hadron and lepton masses and 
the Bd(B s )~meson lifetimes from Ref. [24j . 

All model parameters are fixed by fitting the experimental data in our previous papers (see Refs. [49j - [5ll l54 { ). Their 
numerical values are shown in Eq. (1241) and Table |T| The only new parameter is As which characterizes the size of the 
scalar mesons. We allow this parameter to vary in a relatively large interval, As € [0.8,1.5] GeV. 

In Figs. [3] and Q] we plot our calculated F + (q 2 ) and Fp{q 2 ) form factors in the entire kinematical range 0 < q 2 < q ^ ax . 
Since the behavior of —F-(q 2 ) is very similar to that of F+ (■ q 2 ), we do not display them. One can see that the form 
factors are more sensitive to the choice of As at small q 2 and less so near zero recoil. 

We are going to explore the influence of the intermediate scalar Kq meson on the angular decay distribution of the 
cascade decay B —> Kn + a i + . Therefore, we give the maximum values of the form factors in Tabic Hill and the 
branching ratios in Table HV1 obtained for As = 1.5 GeV. The results for the e mode are almost identical to those of 
the n mode and will not be shown separately. Since the ratio |Vtd|/|Vj s | ss 0.21 is relatively small we do not show 
the branching ratios of the decays with the b — d transition. We compare the obtained results with those from other 
approaches. One can see that our values for the branching ratios are almost half of those from other approaches. 

Let us briefly discuss the impact of the scalar resonance Kq on R —>■ K*(—> K decay. As is well known, the 
narrow AT*(892) vector resonance is described by a Breit-Wigner parametrization and the given cascade B decay can 
be calculated by using the narrow-width approximation. But this is not true in the case of the broad scalar ATp(800) 
meson. There are several parametrizations of the K — n line shapes in the literature; see, for instance, the discussion 
in Ref. For the time being we will use the parametrization accepted in Ref. [r|, the integrated value of which 
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B-S 

F( 0) 

This work 

m 

m 

m 

m 

m 

[55J 

B° d -a+( 980) 

F+( 0) 

0.192 

0.58 

0.56 





B° s -K*+( 800) 

F+( 0) 

0.274 

0.44 

0.53 





B a °-/o(980) 

F+( 0) 

0.254 

0.45 

0.44 

0.19 

0.35 

0.12 

0.40 


F t ( 0) 

0.285 

0.60 

0.58 

0.23 

0.40 

-0.08 


B d ~ Aid 0 (800) 

F+( 0) 

0.306 

0.50 

0.46 






F t ( 0) 

0.306 

0.67 

0.58 





B° d -fo (500) 

F+( 0) 

0.210 








F t ( 0) 

0.203 








TABLE III: The values of the form factors at q 2 = 0 in the covariant quark model (As 


1.5 GeV) and other approaches. 


Decay modes 

Branching fractions 

This work 

(As = 1.5 GeV) 

[37] 

m 

m 

B d —¥ a+ (980)p U/j, 

-)■ dg (980)r P T 

0.52 x 10 -4 

0.11 x 10 -4 

(2.74 ±0.40) x 10 -4 
(1.31 ±0.23) x 10 -4 


1.84 x 10 -4 

1.01 x 10 -4 

B° -»■ K* + ( 800)p"^ 
B° s ->• if*+(800)r P T 

1.23 x 10 -4 

0.25 x 10" 4 

(2.06 ±0.31) x 10" 4 
(1.07 ±0.19) x 10" 4 


1.42 x 10 -4 

0.88 x 10~ 4 

Rd —1 ATg °(800)/i + /r _ 
Rd —>• ATg °(800)r + r“ 

3.47 x 10 -7 

0.61 x 10 -7 

(7.31 ± 1.21) x 10 -7 
(1.33 ±0.36) x 10" 7 



B° s -»■ fo (980)fj, + fi~ 

B° -»■ /o(980)t + t~ 

2.45 x 10" 7 

0.42 x 10 -7 

(5.14 ±0.78) x 10" 7 
(0.74 ±0.17) x 10" 7 

0.95 x 10“ 7 

1.1 x 10~ 7 

5.21 x 10" 7 

0.38 x 10” 7 

Rd -± ATg 0 (800)1?r/ 

5/ /o (980)1?^ 

2.53 x 10 -6 

1.79 x 10~ 6 

(6.30 ±0.97) x 10" 6 
(4.39 ± 0.63) x 10 -6 

0.87 x 10“ 6 



TABLE IV: The branching fractions for the semileptonic and rare B decays into light scalar mesons and lepton pairs. 


in the Af*-resonance region is equal to 


( m K » +8 m ) 2 


( m K » -S m ) 2 


dm 2 Kv \Ls(m 2 K ^.)\ 2 = 0.17, where 6 m = 100 MeV. 


(26) 


Then, we scale the calculated value for the differential decay rate dT(B —> (800)/U + /r ) by this factor and compare 

it with that for B K (892) fi + fi~) decay. We display the behavior of the ratio 


R( 2) = _ 2/3 dT(B —» lT*(892)/i + /r~) _ 

1 2/3 dT(B —»■ K * (892) n + / jl~) + 0.17dT(B —>• K^(800)qi + fi~) 1 ' 

in FigO which may be compared with the finding of Ref. [li|. The integrated ratio (the numerator and denominator 
are integrated separately in the full kinematical region of q 2 ) gives a size for the S- wave pollution to the branching 
ratio of the B —> K*l + l~ decay of about 6%. 































11 


Acknowledgments 

We thank Pietro Santorelli for providing us with the last fit of the parameters in the covariant quark model. We 
would also like to thank Juergen Korner and Valery Lyubovitskij for many useful discussions of S-physics facets 
related to the subject of this paper. We are grateful to Wei Wang and David Straub for pointing out the relevant 
references. 


[1] R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. Ill, 191801 (2013) arXiv:1308.1707 [hep-ex]]. 

[2] R. Aaij et al. [LHCb Collaboration], J. High Energy Phys. 08 (2013) 131 arXiv: 1304.6325 [hep-ex]]. 

[3] T. Hurth and F. Mahmoudi, J. High Energy Phys. 04 (2014) 097 arXiv:1312.5267 [hep-ph]]. 

[4] S. Descotes-Genon, T. Hurth, J. Matias and J. Virto, J. High Energy Phys. 05 (2013) 137 arXiv:1303.5794 [hep-ph]]. 

[5] A. Datta, M. Duraisamy and D. Ghosh, Phys. Rev. D 89, 071501 (2014) arXiv:1310.1937 [hep-ph]]. 

[6] C. Bobeth, G. Hiller, and D. van Dyk, Phys. Rev. D 87, 034016 (2013) arXiv:1212.2321 [hep-ph]]. 

[7] W. Altmannshofer and D. M. Straub, Eur. Phys. J. C 73, 2646 (2013) arXiv:1308.1501 [hep-ph]]. 

[8] W. Altmannshofer and D. M. Straub, State of new physics in b s transitions, arXiv:1411.3161 [hep-ph]]. 

[9] R. Mandal, R. Sinha, and D. Das, Phys. Rev. D 90, 096006 (2014) arXiv: 1409.3088 [hep-ph]]. 

[10] C. D. Lu and W. Wang, Phys. Rev. D 85, 034014 (2012) arXiv:1111.1513 [hep-ph]]. 

[11] M. Doring, U. G. Meissner, and W. Wang, J. High Energy Phys. 10 (2013) 011 arXiv: 1307.0947 [hep-ph]]. 

[12] U. G. Meissner and W. Wang, J. High Energy Phys. 01 (2014) 107 arXiv:1311.5420 [hep-ph]]. 

[13] D. Becirevic and A. Tayduganov, Nucl. Phys. B868, 368 (2013) arXiv: 1207.4004 [hep-ph]]. 

[14] J. Matias, Phys. Rev. D 86, 094024 (2012) arXiv:1209.1525 [hep-ph]]. 

[15] T. Blake, U. Egede, and A. Shires, J. High Energy Phys. 03 (2013) 027 arXiv:1210.5279 [hep-ph]]. 

[16] D. Das, G. Hiller, M. Jung, and A. Shires, J. High Energy Phys. 09 (2014) 109 arXiv:1406.6681 [hep-ph]]. 

[17] A. Ali, P. Ball, L. T. Handoko, and G. Hiller, Phys. Rev. D 61, 074024 (2000) hep-ph/9910221 . 

[18] P. Colangelo, F. De Fazio, P. Santorelli, and E. Scrimieri, Phys. Rev. D 53, 3672 (1996); Phys. Rev. D 57, 3186(E) (1998) 

[19] D. Melikhov, N. Nikitin, and S. Simula, Phys. Rev. D 57, 6814 (1998) 

[20] T. M. Aliev, A. Ozpineci, and M. Savci, Phys. Lett. B 511, 49 (2001) 

[21] C. H. Chen and C. Q. Geng, Nucl. Phys. B636, 338 (2002) 

[22] D. Ebert, R. N. Faustov, and V. O. Galkin, Phys. Rev. D 82, 034032 (2010) arXiv:1006.4231 [hep-ph]]. 

[23] M. A. Ivanov, J. G. Korner, S. G. Kovalenko, and C. D. Roberts, Phys. Rev. D 76, 034018 (2007) 

[24] K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38, 090001 (2014). 

[25] C. Amsler, S. Eidelman, T. Gutsche, C. Hanhart, S. Spanier, and N. A. Tornqvist, Note on Scalar Mesons below 2 GeV 
in Ref. [24^ . 

[26] C. Amsler and N. A. Tornqvist, Phys. Rep. 389, 61 (2004). 

[27] M. Z. Yang, Phys. Rev. D 73, 034027 (2006); Phys. Rev. D 73, 079901(E) (2006) 

[28] C. H. Chen, C. Q. Geng, C. C. Lih, and C. C. Liu, Phys. Rev. D 75, 074010 (2007) 

[29] T. M. Aliev, K. Azizi, and M. Savci, Phys. Rev. D 76, 074017 (2007) arXiv:0710.1508 [hep-ph]]. 

[30] Y. M. Wang, M. J. Aslam, and C. D. Lu, Phys. Rev. D 78, 014006 (2008) arXiv:0804.2204 [hep-ph]]. 

[31] R. H. Li, C. D. Lu, W. Wang, and X. X. Wang, Phys. Rev. D 79, 014013 (2009) arXiv:0811.2648 [hep-ph]]. 

[32] N. Ghahramany and R. Khosravi, Phys. Rev. D 80, 016009 (2009). 

[33] P. Colangelo, F. De Fazio, and W. Wang, Phys. Rev. D 81, 074001 (2010) arXiv: 1002.2880 [hep-ph]]. 

[34] Y. J. Sun, Z. H. Li, and T. Huang, Phys. Rev. D 83, 025024 (2011) arXiv:1011.3901 [hep-ph]]. 

[35] H. Y. Han, X. G. Wu, H. B. Fu, Q. L. Zhang, and T. Zhong, Eur. Phys. J. A 49, 78 (2013) arXiv:1301.3978 [hep-ph]]. 

[36] Z. G. Wang, Eur. Phys. J. C 75, 50 (2015) arXiv: 1409.6449 [hep-ph]]. 

[37] Z. G. Wang, Semi-leptonic B —>- S decays in the standard model and in the universal extra dimension model, 
arXiv:1411.7961 [hep-ph]]. 

[38] A. J. Buras and M. Munz, Phys. Rev. D 52, 186 (1995) hep-ph/9501281 . 

[39] G. Buchalla, A. J. Buras, and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996) hep-ph/9512380 . 

[40] A. Ali, T. Mannel, and T. Morozumi, Phys. Lett. B 273, 505 (1991). 

[41] A. Khodjamirian, T. Mannel, A. A. Pivovarov, and Y.-M. Wang, J. High Energy Phys. 09, (2010) 089 arXiv: 1006.4945 
[hep-ph]]. 

[42] A. Faessler, T. Gutsche, M. A. Ivanov, .J. G. Korner, and V. E. Lyubovitskij, 

Eur. Phys. J. direct C 4, 18 (2002) hep-ph/0205287 . 

[43] F. Kruger and J. Matias, Phys. Rev. D 71, 094009 (2005) hep-ph/0502060 . 

[44] T. Branz, A. Faessler, T. Gutsche, M. A. Ivanov, J. G. Korner, and V. E. Lyubovitskij, Phys. Rev. D 81, 034010 (2010) 
arXiv:0912.3710 [hep-ph]]. 

[45] S. Weinberg, Phys. Rev. 130, 776 (1963). 

[46] A. Salam, Nuovo Cimento 25, 224 (1962). 

[47] K. Hayashi, M. Hirayama, T. Muta, N. Seto, and T. Shirafuji, Fortsch. Phys. 15, 625 (1967). 

[48] G. V. Efimov and M. A. Ivanov, The Quark Confinement Model of Hadrons, (CRC Press, Boca Raton, 1993). 



12 


[49] M. A. Ivanov, J. G. Korner, S. G. Kovalenko, P. Santorelli, and G. G. Saidullaeva, Pliys. Rev. D 85, 034004 (2012) 
arXiv:1112.3536 [hep-ph]]. 

[50] S. Dubnicka, A. Z. Dubnickova, M. A. Ivanov, and A. Liptaj, Phys. Rev. D 87, 074201 (2013) arXiv:1301.0738 [hep-ph]]. 

[51] T. Gutsche, M. A. Ivanov, J. G. Korner, V. E. Lyubovitskij, and P. Santorelli, Phys. Rev. D 86, 074013 (2012) 
arXiv: 1207.7052 [hep-ph]]; Phys. Rev. D 87, 074031 (2013) arXiv:1301.3737 [hep-ph]]; Phys. Rev. D 88, no. 11, 114018 
(2013) arXiv:1309.7879 [hep-ph]]; Pliys. Rev. D 90, no. 11, 114033 (2014) arXiv:1410.6043 [hep-ph]]. 

[52] S. Dubnicka, A. Z. Dubnickova, M. A. Ivanov, and J. G. Korner, Phys. Rev. D 81, 114007 (2010) arXiv: 1004.1291 [hep- 
ph]]; S. Dubnicka, A. Z. Dubnickova, M. A. Ivanov, J. G. Koerner, P. Santorelli, and G. G. Saidullaeva, Phys. Rev. D 84, 
014006 (2011) arXiv: 1104.3974 [hep-ph]]. 

[53] C. Y. Cheung and C. W. Hwang, J. High Energy Phys. 04, (2014) 177 arXiv:1401.3917 [hep-ph]]. 

[54] G. Ganbold, T. Gutsche, M. A. Ivanov, and V. E. Lyubovitskij, On the meson mass spectrum in the covariant confined 
quark model, arXiv:1410.3741 [hep-ph]]. 

[55] B. El-Bennich, O. Leitner, J.-P. Dedonder, and B. Loiseau, Phys. Rev. D 79, 076004 (2009) arXiv:0810.5771 [hep-ph]]. 


13 







FIG. 3: The F+ ( q 2 ) form factors for the b — u, b — d, and b — s transitions. The upper and lower edges correspond to the values 
A s = 1.5 GeV and As = 0.8 GeV, respectively. 
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FIG. 4: The Fr{q 2 ) form factors for the b — u, b — d, and b — s transitions. The upper and lower edges correspond to the values 
As = 1.5 GeV and As = 0.8 GeV, respectively. 



FIG. 5: The ratio of the differential decay rate dF(B —> K*( 892)(—> K 0 n + )fj, + fj, )) to the full differential decay rate dF(£> — ► 
A'*(892)(-> A+7r-)/i+ M -) + dF{B ->• K+n-)fj+fj~). 






















